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We investigate static axially symmetric monopole and black hole solutions with magnetic charge 
n > 1 in Einstein- Yang-Mills-Higgs theory. 

For vanishing and small Higgs selfcoupling, multimonopole solutions are gravitationally bound. 
Their mass per unit charge is lower than the mass of the n = 1 monopole. For large Higgs selfcoupling 
only a repulsive phase exists. 

The static axially symmetric hairy black hole solutions possess a deformed horizon with constant 
surface gravity. We consider their properties in the isolated horizon framework, interpreting them 
as bound states of monopoles and black holes. Representing counterexamples to the "no-hair" 
conjecture, these black holes are neither uniquely characterized by their horizon area and horizon 
charge. 
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I. INTRODUCTION 



Magnetic monopoles |l[ , multimonopoles and nionopole-antimonopole pair solutions are globally regular 
solutions of SU(2) Yang-Mills-Higgs (YMH) theory, with Higgs field in the triplet representation. Since their magnetic 
charge is proportional to their topological charge, the monopoles and multimonopoles reside in topologically non-trivial 
sectors of the theory, whereas the monopole-antimonopolc pair solutions are topologically trivial. 

In the Bogomol'nyi-Prasad-Sommerfield (BPS) limit where the strength of the Higgs self-interaction potential 
vanishes, the mass of the monopole and multinionopole solutions saturates its lower bound, the Bogomol'nyi bound. 
In particular, the mass per unit charge of an n > 1 monopole is precisely equal to that of the n — 1 monopole. The 
massless Higgs field mediates a long range attractive force which exactly cancels the long range repulsive magnetic 
force of the U{1) field 

For finite Higgs selfcoupling, however, the Higgs field is massive and therefore decays exponentially. Consequently 
the long range magnetic field dominates at large distances, leading to the repulsion of like monopoles [|2|. In particular, 
as verified numerically for n = 2 and n = 3 monopoles the mass per unit charge of an n > 1 monopole is higher 
than the mass of the n = 1 monopole. Thus for finite Higgs selfcoupling there is only a repulsive phase between like 
monopoles. 

Let us now consider the effect of gravity on the monopole and multinionopole solutions. When gravity is coupled 
to YMH theory, a branch of gravitating monopole solutions emerges from the flat space monopole solution . With 
increasing gravitational strength, the mass of the gravitating monopole solutions decreases monotonically. The branch 
of monopole solutions extends up to some maximal value of the gravitational strength, beyond which the size of the 
soliton core would be smaller than the Schwarzschild radius of the solution . The same holds true for static axially 
symmetric gravitating multimonopole solutions ]T^ . All these solutions are asymptotically flat. 

The inclusion of gravity allows for an attractive phase of like monopoles not present in flat space . There arises 
a region of parameter space, where the mass per unit charge of the gravitating multimonopole solutions is lower than 
the mass of the gravitating n = 1 monopole. Here the multimonopole solutions are gravitationally bound. 

To every regular monopole solution there exists a corresponding family of black hole solutions with regular event 
horizon and horizon radius < xa < 2:A,max fl^- Likewise, to every regular axially symmetric multimonopole solution 
there exists a corresponding family of black hole solutions with regular event horizon p4| . Outside their event horizon 
these black hole solutions possess non-trivial non-abelian fields. Therefore they represent counterexamples to the 
"no-hair" conjecture. The axially symmetric black hole solutions additionally show, that static black hole solutions 
need not be spherically symmetric, i. e. Israel's theorem cannot be generalized to non-abelian theories either [p"5|-p7|. 

Considering the non-abelian black hole solutions in the isolated horizon framework |18-20j, they can be interpreted 



as bound states of monopoles and black holes pO[ . In particular, the isolated horizon framework yields an intriguing 
relation for the mass of hairy black hole solutions, representing it as the sum of the monopole mass and the horizon 
mass of the black hole solutions |^ . Having shown previously, that this relation is also valid for black holes inbetween 
monopole-antimonopolc pair solutions we here verify this relation for the magnetically charged hairy black hole 
solutions. The isolated horizon formalism has furthermore led to new conjectures for black holes. In particular, a 
"quasilocal uniqueness conjecture" has been proposed, stating that static black holes are uniquely characterized by 
their horizon area and horizon charge(s) We investigate the validity of this conjecture for the EYMH black hole 
solutions. 

This paper presents a detailed account of the static axially symmetric multimonopole and black hole solutions, 
reported in p4[| . In section II we present the action, the axially symmetric ansatz in isotropic spherical coordinates 
and the boundary conditions. In section III we recall the spherically symmetric solutions, presenting them in isotropic 
coordinates. In section IV we discuss the properties of the axially symmetric regular multimonopole solutions, and in 
section V those of the black hole solutions. We present our conclusions in section V. In Appendix A some details of 
the quantities -F)ii/, -0^$ and T^j^ are shown, and in Appendix B the numerical technique is briefly described. 

II. EINSTEIN- YANG-MILLS-HIGGS EQUATIONS OF MOTION 
A. Einstein- Yang-Mills-Higgs action 

We consider the SU(2) Einstein- Yang-Mills-Higgs action 



^ " / ( 16^G " iTr(F^,F^'^) - ^Tr(i?^$i5^<f>) - |ATr($2 _ rj^f ) (1) 



with field strength tensor 
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of the gauge field 



and with covariant derivative 



F^^ = df,A„ - d^Af, + ie [A^, A^\ , (2) 



A, - ^t^A; , (3) 



D^<^> = d^^ + ie[A^,<^>] , (4) 
of the Higgs field in the adjoint representation 

$ = . (5) 

Here g denotes the determinant of the metric. The constants in the action represent Newton's constant G, the Yang- 
Mills coupling constant e, the Higgs self-coupling constant A and the vacuum expectation value of the Higgs field 
V- 

Variation of the action (|l|) with respect to the metric g'^'^ leads to the Einstein equations 

Gpi/ — -R/ji/ ~ -^gfivR = SirGTf^^ (6) 

with stress-energy tensor 



Tfii, — gi_ii/L}^i — 2 



= Tr(iz?^$i?,$ - ig^,i?„$i?"$) + 2Tr(F^,F,/35"'^ " 75^.^^0/3^"'^) - \g^M{<^^ - if f . (7) 
z 4 4 o 

Variation with respect to the gauge field and the Higgs field $ leads to the matter field equations, 

^ .D^i^Fn - 7*e[$, D^^ = , (8) 



'A 

^D^X^D^'^) + A($2 _ ^f)^ ^ , (9) 



respectively. 



B. Static axially symmetric ansatz 



Instead of the Schwarzschild-like coordinates, used for the spherically symmetric EYM and EYMH solutions 
p2| , p3|p3| (see section HI), we adopt isotropic coordinates as in | |2^ , p^ , p5| , pj| . p6p]]Jl^ , to construct static axially 
symmetric solutions. In terms of the spherical coordinates r, 9 and the isotropic metric reads 

,9 ,. , 9 , 9 JTir^ Ir^ sin^ 6* , 9 , 
ds^ = -fdf + jdr^ + + J dip^ , (10) 

where the metric functions /, m and I are only functions of the coordinates r and 6. The z-axis {9 = 0) represents 
the symmetry axis. Regularity on this axis requires 

TO|e=o = l\0=o ■ (11) 
We take a purely magnetic gauge field, ^0 — 0, and choose for the gauge field the ansatz I^ j2^ , p^ , |l6| , p5| , p7 14 



A^dx^" = — [t; {Hidr + (1 - iJa) rdO) - n {t'^H^ +T]}{l-Hi))r sin Odp^] . (12) 

Here the symbols r", Tg and t"^ denote the dot products of the cartesian vector of Pauli matrices, r = {Tx,Ty,Tz), 
with the spatial unit vectors 
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e " = (sin 9 cos rup, sin 9 sin mp, cos 9) 
eg = (cos 9 cos n(/?, cos 6* sin nip, ~ sin ( 



(— sinniyj, cosn(/3, 0) , 



(13) 



respectively. Since the fields wind n times around, while the azimuthal angle (p covers the full trigonometric circle 
once, we refer to the integer n as the winding number of the solutions. For the Higgs field the corresponding ansatz 



IS 



il 



(14) 



The four gauge field functions Hi and the two Higgs field functions <i>i depend only on the coordinates r and 9. For 
71 = 1 and Hi — = ^2 = 0, H2 = H4 = K{r) and $1 = H{r), the spherically symmetric solutions are obtained in 
isotropic coordinates. 

The ansatz (p^)-([l^) is axially symmetric in the sense, that a rotation around the z-axis can be compensated by a 
gauge rotation. The ansatz is form-invariant under the abelian gauge transformation |0,^,^,Eol 



(15) 



C/ = exp(^-T;T(r,( 

The functions Hi and H2 transform inhomogeneously under this gauge transformation, 

Hi^ Hi- rdrT , 
H2^H2+ deV , 



(16) 



like a 2-dimensional gauge field. The functions H3 and H4 combine to form a scalar doublet, (iJa+cotfl, H4). Likewise, 
the Higgs field functions form a scalar doublet ($1, $2)- 

We fix the gauge by choosing the gauge condition as previously ||2^,^,^,^,^ 25 l^Jl^ . In terms of the functions 
Hi it reads 



rdrHi - dgH2 = . 

With the ansatz (p^-(|T^) and the gauge condition (|l^) we then obtain the set of EYMH field equations. 
The energy density of the matter fields e = —Tq = —Lm is given by 



(17) 



^0 ^ Try^J 



2r2 



- {{rdr<^i + Hi^2f + {rdr'^'2 - Hi^if + (ae$i - i/2$2)^ + {de^2 + i?2$i) 



- — (iJ4$l + (i^3+COt0)$2)' 



p 



- {rdrH2 + deHiY 



{rdrHs - HiHif + {rdrHi + Hi {H3 + coi9)f 



+ {dgH3 - 1 + C0i9Hi + HiHif + {dgH4 + cot9 {H4 - H2) - i?2i?3)'] } 
As seen from Eq. (p^), regularity on the z-axis requires 

H2\e=Q ~ -ff4|e=o • 



(18) 



(19) 



C. Boundary conditions 



To obtain asymptotically flat solutions with the proper symmetries, which are either globally regular or possess 
a regular event horizon, we must impose appropriate boundary conditions p^ , p^p5 T^,]!^. Here we are looking for 
solutions with parity reflection symmetry. Therefore we need to consider the solutions only in the region Q < 9 < ^, 
imposing boundary conditions along the p- and z-axis (i. e. for 9 — ^ and 9 = Q). 

Boundary conditions at infinity 
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Asymptotic flatness imposes for the metric functions of the regular and black hole solutions at infinity (r = oo) the 
boundary conditions 

f\r—oo — ^|r— oc — ^|r— oo — 1 • (^0) 

For the Higgs field functions we require 

*l|r=oo = l, $2|r=oo-0, (21) 

thus the modulus of the Higgs field assumes the vacuum expectation value rj. For magnetically charged solutions, the 
gauge field functions Hi satisfy 

i?l|r=oo — H2\r=oo = -ff3|r=oo = ^^4|r=oo — . (22) 

Boundary conditions at the origin 

Requiring the solutions to be regular at the origin (r = 0) leads to the boundary conditions for the metric functions 

drf\r=0 = 9,.m|,.=o = drl\r=0 = . (23) 

The Higgs field functions satisfy 

$l|r=0 = *2|r=0 = , (24) 

and the gauge field functions Hi satisfy 

i?l|r=0 — Hslr^Q = , H2\r=0 — i?4|r=0 1 ■ (25) 

Boundary conditions at the horizon 

The event horizon of static black hole solutions is characterized by gu = — / = 0. In isotropic coordinates grr is 
finite at the horizon. We now impose that the horizon of the black hole solutions resides at a surface of constant r, 

Requiring the horizon to be regular, we obtain the boundary conditions at the horizon r — vn. The metric functions 
must satisfy 

/|r=rH = m\r=rii = l\r=ni = . (26) 

The boundary conditions for the gauge field functions and the Higgs field functions can be deduced from the equations 
(^ and (^), respectively. 

Fr9\r=r^ = ^ (rS^i/a + dgHi)\r=r^ = , 

f^r^|r=rH = ^ {rOrHs - HiHi)\r=r^ = , {rdrH^ - Hi{H:i + COt0))|,=,„ = , 
Dr<^\r=r^=Q ^ {vdr^l + Hi^2)\r=r^ ^ ^ . {rdr^2 - Hi^ i)\r=r^ = Q . (27) 

The equations of motion yield only three boundary conditions for the four gauge field functions Hi\ one gauge field 
boundary condition is left indeterminate. However, for the black hole solutions the gauge condition (^7|) still allows 
non-trivial gauge transformations satisfying 

r'^dlT + rdrT + 5|r = . (28) 

To fix the gauge, we have chosen the gauge condition ||l^Jl7| , p!^ 

{deHi)\r^r^ = 0, (29) 

which implies -ffi|r=rH = Oi if we take into account the boundary conditions on the axes, -ffi|e=o,f = Oi (see next 
paragraph). As a consequence the boundary conditions Eqs. ( |27| ) reduce to 

drH2\r=r^ = 0, drH:i\r=r^ = 0, drHi\r=r^ = 0, c',*!!,^,^ - 0, dr^2\r=r^ = 0. (30) 

Boundary conditions along the axes 

The boundary conditions along the p- and z-axis (6 ~ ^ and 6 — Q) are determined by the symmetries. The metric 
functions satisfy along the axes 
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d9f\e=o = dem\g=o = dglle^o = , 

def\e=^ = dgm\e^^ = del\e=^ ^ . (31) 
Likewise the Higgs field functions satisfy 

9e$i|e=o = 0, $2|0=o = O, 

9e*i|e=f=0, ^2\e=^=0, (32) 
along the axes. For the gauge field functions Hi symmetry considerations lead to the boundary conditions 

Hi\o=o = i?3|e=o = , d0H2\e=o — 9ei?4|e=o = , 

Hi\e=i = H3\9=i = , dgH2\e=i = d9Hi\e=^ = , (33) 

along the axes. In addition, regularity on the z-axis requires condition (^l]) for the metric functions to be satisfied, 
and condition ( [l9| ) for the gauge field functions. 

D. Mass, temperature and entropy 

Let us introduce the dimensionless coordinate x, 

X — rjer . (34) 
The equations then depend only on two dimensionless coupling constants, a and /3, 

^ ATTGif , = 4 ■ (35) 

The mass M of the solutions can be obtained directly from the total energy-momentum "tensor" t^'^ of matter and 
gravitation 

M = j T°°d^r . (36) 

It is related to the dimensionless mass n/a^, via 

^./a^ = , (37) 



where /i is determined by the derivative of the metric function / at infinity |g^,|l6|,^ 17 

/i = i hm x'^dxf ■ (38) 

Z X — >oo 

The surface gravity of static black hole solutions is given by ||3^Jl^ 

4 = "(l/4)ffV(5.5«)(a,g„) . (39) 

To evaluate Hsg, we need to consider the metric functions at the horizon. Expanding the equations in the vicinity of 
the horizon in powers of the dimensionless coordinate {x — x-ii)/x-ii, we observe, that the metric functions are quadratic 
in a; — XH, 

= MO) ('-^y (i - ^) + o (^-^y , (40) 

\ xn J \ xa J \ xn J 

m(x,e) = m^ie) ( ^—^y ( 1 - + O ( ^—^Y , (41) 

with l{x, 9) like m{x, 9), Eq. (p]). We then obtain the dimensionless surface gravity k — Ksg/erj 
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f2{0) 



(42) 



The zeroth law of black hole physics states, that the surface gravity Kgg is constant at the horizon of a black hole 
[Q. To show that the zeroth law holds for the hairy black hole solutions we employ the expansion of the metric 
functions (^o|)-(^l|) in the r9 component of the Einstein equations at the horizon [0. This yields the crucial relation 
between the expansion coefficients /2(^) and m2{0), 







dem2 _ ^ ^9/2 

TO2 /2 



The temperature is proportional to the surface gravity Kg 

by 



T = k/{2tt) . 

The dimensionless area A of the event horizon of the black hole solutions 



(43) 

in particular the dimensionless temperature is given 

(44) 



A = 271 

is proportional to the dimensionless entropy S [p! 



,a ■ nVh^ 2 
dO sm fc* — xfj 

72 



yielding the product 



TS = 



d6 sineWl 



(45) 



(46) 



(47) 



Having defined temperature and entropy, we now derive a second expression for the mass of the black hole solutions 
52| . As in EYM theory , the equations of motion yield in EYMH theory the relation 



-g - T/) 



Integrating both sides over r, 9 and ip from the horizon to infinity, we obtain 

/>2tt />n noo 



1 r 

— / dOsind 
4G Jo 



nZTV nTV nOQ 

/ / / d^d0drV^(2ro*'-T/) =Mo 

Jo Jo JrH 



(48) 



(49) 



Changing to dimensionless coordinates, we express the l.h.s. with help of the dimensionless mass /x and the product 
of temperature and entropy, TS*, obtaining 



H = Ho + 2TS , 



(50) 



with iioja^ 



Afo, in agreement with the general mass formula for static black hole solutions Isa]. For regular 



solutions one simply obtains 



M 



(2To° - T/) 



-gdrdOdip 



We finally consider the Kretschmann scalar K, 



(51) 



(52) 



In order to show that K is finite at the horizon it is suffient to point out, that the expansion of the metric functions 
(^)-(P) is the same as for the EYM system considered in |17|. Consequently, it follows from the calculations of 
ref. ]17[ that the Kretschmann scalar is finite at the event horizon, if the condition (^) holds, i. e. if the temperature 
is constant at the event horizon. 
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E. Horizon mass and horizon charge 



In the isolated horizon framework ||T^ an intriguing relation between the ADM mass n/a^ of a black hole with area 
A and area parameter xa , 



a;A - ^/A/A^T , (53) 
and the ADM mass /iiog/a^ of the corresponding regular solution holds |p^ , 

= Mrog + MA , (54) 

where the horizon mass [it^lo? is defined by 



^0 



(55) 



with the dimensionless surface gravity k in the integrand. 

The isolated horizon formalism then suggests to interpret a hairy black hole as a bound state of a regular solution 
and a Schwarzschild black hole | pO| , 

H = Mrcg + /is + /-ibind , (56) 

where ^is/oP' — x/^/2a^ is the ADM mass of the Schwarzschild black hole with horizon radius xa, and /ibind/a'^ 
represents the binding energy of the system, 

Mbind = ^J■A- ^J■s ■ (57) 



Another crucial quantity in the isolated horizon formalism is the non-abelian magnetic charge of the horizon |18| , |19[ 
defined via the surface integral over the horizon 



The non-abelian electric horizon charge is defined analogously with the dual field strength tensor |1^J19[]. These 
horizon charges are an important ingredient in a new "quasilocal uniqueness conjecture" put forward in ref. p9| , 
which states that static black holes are uniquely determined by their horizon area and their horizon electric and 
magnetic charges. 

III. SPHERICALLY SYMMETRIC SOLUTIONS 

Previously, spherically symmetric EYMH solutions were obtained numerically in Schwarzschild-like coordinates 
p3[ . Since we construct the axially symmetric solutions in isotropic coordinates, we here discuss the coordinate 
transformation between isotropic and Schwarzschild-like coordinates for spherically symmetric solutions. 

We briefly recall the dependence of the monopole solutions on the parameters of the theory, a and /?. In particular, 
we compare the way the limiting RN solution is reached in Schwarzschild-like and in isotropic coordinates. 

We then turn to the spherically symmetric black hole solutions, putting particular emphasis on the relations obtained 
in the isolated horizon formalism. We demonstrate the relation between black hole mass, soliton mass and horizon 



mass, Eq. (p4), we evaluate the binding energy, Eq. (56), and we discuss the "quasilocal uniqueness conjecture" [Il9l. 



A. Coordinate transformations 

By requiring I — m and the metric functions / and m to be only functions of the coordinate r, the axially symmetric 
isotropic metric (^) reduces to the spherically symmetric isotropic metric 

ds^ = -fdt^ + J {dr^ + {dd^ + sin^ ed<fi^)) . (59) 
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In Schwarzschild-like coordinates the metric reads |lc 



1 



where the metric functions A and N are only functions of the radial coordinate f, and 

iV(f).l-^. 



(60) 



(61) 



The spherically symmetric horizon resides at radial coordinate fu , and the (dimensionfuU) area of the horizon is 

A = Anr^ , (62) 

hence the area parameter ta is just the Schwarzschild radius ffj. 

Comparison of the metric in Schwarzschild-like coordinates ( |6C| ) with the isotropic metric ( |59| ) yields for the coor- 
dinate transformation 



1 df 



dr 



(63) 



The function N{f) (or equivalently the mass function M{f)) is only known numerically. Therefore the coordinate 
function r{f) can only be obtained numerically from Eq. (p3[). 



B. Monopole solutions 



Let us briefly recall the dependence of the gravitating magnetic monopole solutions on the parameters a and (3. 
When a is increased from zero, while (3 is kept fixed, a gravitating monopole branch emerges smoothly from the 
corresponding flat space monopole solution. This fvmdamental monopole branch extends up to a maximal value amax 
of the coupling constant a. Beyond this value no monopole solutions exist. In the BPS limit, /3 = 0, and for small 
values of /?, the fundamental monopole branch bends backwards at Ofmax, until a critical coupling constant acr is 
reached. At the critical value the fundamental monopole branch reaches a limiting solution and bifurcates with 
the branch of extremal RN solutions of unit magnetic charge ||l^ . For larger values of /3 the maximal value of a and 
the critical value of a coincide, amax = acr [p^ 33 3^ . 

Let us first recall, how the critical value of a is approached in Schwarzschild-like coordinates. Along the fundamental 
branch, the metric function N{x) of the monopole solutions develops a minimum, which decreases monotonically. In 
the limit a — > acr, the minimum approaches zero at Xcr = ckcr- The limiting metric function then consists of an inner 
part, X < Xcr, and an outer part, x > icr- For x > Xct, the limiting metric function corresponds to the metric function 
N{x) of the extremal RN black hole for acr with unit magnetic charge. 

Likewise the gauge and Higgs field functions approach limiting functions, when a — » acr- For x > icr they 
also correspond to those of the extremal RN black hole for acr with unit magnetic charge. The limit a —^ acr is 
demonstrated in Figs, la-c for the monopole solutions in the BPS limit for the metric function N{x), the gauge field 
function K{x) and the Higgs field function H{x). 

We recall, that a RN solution with magnetic charge P has metric functions 



N{i) = 1 - ^ + . Mi) = 1 • 



(64) 



In the embedded RN solution with unit magnetic charge the gauge field function K{x) and the Higgs field function 
H{x) are constant. 



K{3:) = , H{3:) = 1 



(65) 



implying a Coulomb-type decay for the magnetic field and a trivial Higgs field, assuming its vacuum expectation 
value. In particular, an extremal RN solution has metric function N{x), 



N{i) 



aP 

X 



(66) 
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and ADM mass y^jo? ^ where 



M = a|P| . (67) 

Let us now turn to the monopole solutions in isotropic coordinates. In Figs. 2a-c we demonstrate the dependence of 
the gravitating magnetic monopole solutions on the parameter a along the fundamental branch in isotropic coordinates. 
For comparison the same set of parameter values as in Figs, la-c is chosen. In the limit a acr the fundamental 
monopole branch bifurcates with the branch of extremal RN solutions of unit magnetic charge. In particular, as the 
critical a is approached, the metric function f{x) develops a zero at the origin, which corresponds to the horizon of 
an extremal RN solution in isotropic coordinates. 

For a RN solution with charge P and horizon radius the isotropic coordinate x is related to the Schwarzschild-like 
coordinate x by 



JxP' — lax + cP-P"^ + X — u 

n (68) 



in the non-extremal case a\P\ < /i = (a;^ + a^P^), and by 

2xh 



■a\P\ (69) 



in the extremal case a\P\ = ii ~ xn. 

The metric functions f{x) and m{x) of a non-extremal RN black hole are given by 

fix) - ^ ^ xn' 



mix) = (^)\l-^ni + ^r, (70) 
\ x / xh Xfi 

respectively. The horizon radius xn in isotropic coordinates is related to the horizon radius in Schwarzschild-like 
coordinates xa by 

xl-a'P' ^71^ 

= 7-2 ■ (71) 



4xi 



Its ADM mass /i/a^ is obtained from 



An extremal RN solution has horizon radius xn = and metric functions 

2 

X + a\P 



/(^)-( , m(x) = l. (73) 



C. Black hole solutions 



Let us now turn to the black hole solutions of the SU(2) EYMH system. We here limit our discussion to the BPS 
case, (3 = 0. 

Hairy black hole solutions exist in a limited domain of the xa — a-plane |]l3| . For fixed a < amax, hairy black hole 
solutions emerge from the monopole solution in the limit xa 0. They persist up to a maximal value of the horizon 
radius a;A,max, which limits the domain of existence of hairy black holes. The domain of existence is shown in Fig. 3, 
where the maximal value of the horizon radius XA.max is shown as a function of a. 

The domain of existence of hairy black hole solutions consists of two regions with distinct critical behaviour. These 
two regions are separated by the particular value of a, a = ^3/2 [|l^. For a < a < ttmax, the hairy black hole 
solutions bifurcate at a critical value XA,cr with an extremal RN solution with unit charge. In contrast, for < a < ci, 
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the hairy black hole solutions bifurcate at a critical value xa.ci: with a non-extremal RN solution with unit charge. In 
particular, for small values of a, the branch of hairy black hole solutions extends backwards from the maximal value 
a^A.max to the critical value XA,cr, whereas for larger values of a < a, both values coincide, XA.max — a^A.cr- 

The dependence of the mass fi/a^ of the hairy black hole solutions on the area parameter is demonstrated in Fig. 4 
for a = 0.5 and 1. For a = 0.5 < d, the mass of the black hole solutions increases monotonically with increasing xa, 
until it reaches its maximal value at XA,ma.x- Bending backwards the mass then decreases again, until the bifurcation 
point at XA,ci- is reached, where it coincides with the mass of a non-extremal RN black hole with unit charge. Note, 
that the mass of the hairy black hole solutions exceeds the mass of the RN solutions in a small region close to XA.max- 

In contrast for a = 1 > d, the limiting solution reached corresponds to an extremal RN solution with unit charge 
(for X > XA.cr)- Since for a > a the maximal horizon radius a;A,inax of the hairy black holes is smaller than the 
horizon radius of the corresponding extremal RN solution, a gap between the branch of hairy black holes and the RN 
branch arises. This gap is seen in Fig. 4 for a = 1 for the mass of the black holes. 

In Fig. 5 we exhibit the surface gravity k as function of the area parameter xa for the same set of black hole 
solutions, obtained for a = 0.5 and 1. The surface gravity of the hairy black hole solutions decreases monotonically 
along the black hole branches. In contrast, the surface gravity of the corresponding RN branches increases for small 
horizon radius xa. For xa a^A,cr, the surface gravity of the hairy black hole solutions reaches the surface gravity 
of the corresponding limiting RN solutions. In particular, for a ~ 0.5 the hairy black hole branch and the RN branch 
bifurcate, and the limiting value reached corresponds to the value of a non-extremal RN solution. In contrast for 
a = 1 the surface gravity of the hairy black hole solutions reaches the value of an extremal RN solution, namely zero, 
even though the hairy black hole branch is separated from the RN branch by a gap of the horizon radius. 

Let us now turn to the predictions of the isolated horizon formalism, which have not been considered before for 
EYMH black hole solutions. We first consider relation ( p^ ) between black hole mass, soliton mass and horizon mass for 
the fundamental monopole and black hole solutions. In Fig. 6 we exhibit the horizon mass [iaIcP' obtained from js^). 
Adding the mass of the corresponding soliton solution (/Ltreg/a^(a = 0.5) = 0.963461 and ^rcg/a^(a = 1-0) = 0.855254, 
respectively,) precisely gives the masses of the non-abelian black hole solutions shown in Fig. 4. 

Interpreting the hairy black holes as a bound states of regular solutions and Schwarzschild black holes , we can 
identify the binding energy of these systems /ibind/o^^, according to Eq. (|5^). In Fig. 7 wc present the binding energy 
of the hairy black hole solutions for a — 0.5 and a = 1. For comparison, we also show the binding energy MbindV'^^ 
of the RN solutions, which we define analogously to Eq. (|5^ ) via 

Mbind = MRN - Mrcg - MS , (74) 

where /iRN/a^ is the mass of the RN solution with unit charge, /Zrog/o;^ is the mass of the monopole, and o? is the 
mass of the Schwarzschild black hole with the same horizon area as the RN black hole. Note that, when the solutions 
coexist, the binding energy of the hairy black hole solution is smaller than the binding energy of the RN solution 
(except for a small region close to XA.max (see also Fig. 4)). This indicates stability of these hairy black hole solutions 
(except close to a;A,max)- 

Let us finally consider the "quasilocal uniqueness conjecture" of ref. which states, that static black holes are 
uniquely specified by their horizon area and horizon charges. The solutions considered here carry no horizon electric 
charge, thus they should be uniquely characterized by their horizon area and horizon magnetic charge. In Fig. 8 we 
exhibit the non-abclian horizon magnetic charge for the hairy black hole solutions with a = 0.5 and 1. The horizon 
magnetic charge increases monotonically along the branches of hairy black hole solutions. As expected, for a — 0.5 
the value of the RN solution with unit charge is reached, when xa xa.ct- In contrast, for a = 1 a value smaller 
than one is reached, when xa xa.ci- In this case again a gap occurs between the hairy black hole branch and the 
RN branch. 

Concerning the "quasilocal uniqueness conjecture" we conclude, that allowing only for integer values of the magnetic 
charge as required for the non-abelian magnetic charge of embedded RN solutions, the spherically symmetric black 
hole solutions are uniquely characterized by their area parameter and their horizon magnetic charge. 

IV. AXIALLY SYMMETRIC SOLUTIONS 

We here give a detailed discussion of the properties of the regular axially symmetric multimonopole solutions. 
We first investigate the dependence of the monopole solutions on the parameter a for fixed (3. In particular, we 
demonstrate the convergence of the multimonopole solutions with magnetic charge n to limiting extremal RN solutions 
with magnetic charge n, when a —^ ctcr. Focussing on values of a close to ctcr, we introduce auxiliary Schwarzschild-like 
coordinates to gain better understanding of the limiting solutions obtained in isotropic coordinates. 
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We then show, that the inclusion of gravity allows for an attractive phase of like monopoles not present in flat space 
1^^ . There arises a region of parameter space, where the mass per unit charge of the gravitating multimonopole solu- 
tions is lower than the mass of the gravitating n — 1 monopole, hence the multimonopole solutions are gravitationally 
bound. 

The numerical technique is briefly described in Appendix B. 



A. Fundamental multimonopole branch 



Let us first consider the dependence of the gravitating axially symmetric multimonopole solutions with magnetic 
charge n on the parameters a and /3. Analogously to the monopole solutions, when a is increased from zero, while f3 is 
kept fixed, a branch of gravitating multimonopole solutions with charge n emerges smoothly from the corresponding 
flat space multimonopole solution. The multimonopole branch extends up to a maximal value aniax('^) of the coupling 
constant a, beyond which no axially symmetric multimonopole solutions with charge n exist. At the maximal value 
ctmaxin), which coincides with the critical value acj:{n) [|35| , the multimonopole branch reaches a limiting solution 
and bifurcates with the branch of extremal RN solutions with magnetic charge n jl^,^ . The metric functions of this 
embedded RN solution are given by Eqs. (73) with I = m, and the gauge field and Higgs field functions are constant, 
Hi = 0,i = l, ...,4, $1 = 1, and $2 = 0. 

We now illustrate the dependence of the multimonopole solutions on a, and in particular the convergence of the non- 
abelian solutions to the corresponding RN solution in the limit a acr, for the special case of n = 2 multimonopole 
solutions in the BPS limit. Numerically we find in this case acr ~ Q^max ~ 1.5 |p^ , ^ . The static axially symmetric 
solutions depend on two variables, the radial coordinate x and the angle 9. In the following we present the functions 
in two-dimensional plots, exhibiting the x-dependence for three fixed angles, 9 = 0, 9 = tt/A and 9 — Tr/2. 

Let us first discuss the dependence of the metric functions on the parameter a. In Figs. 9a-c the metric functions /, 
I and TO are shown, respectively, for a = 1, 1.2, 1.4 and 1.499. The function / increases monotonically with increasing 
X for all values of a. Its value at the origin /(O) decreases with increasing a and tends to zero as a approaches the 
critical value a^- The functions ra and I also increase monotonically with increasing x. As a tends to its critical 
value, these functions approach the value one on an increasing interval. However, at the origin the functions to and I 
assume a value different from one. Thus, convergence to the limiting RN solution is pointwise. 

In Figs. 9d-g we show the gauge fields functions H1-H4, respectively, for the same values of a. The function Hi 
possesses a maximum, whose position decreases with increasing a and tends to zero, when a approaches its critical 
value. In contrast, the height of the maximum depends only weakly on a. The function H2 decreases monotonically 
with increasing x on the z-axis {9 = 0) and possesses a minimum on the p-axis (9 = 7r/2). The location of the 
minimum decreases with increasing a and tends to zero when a approaches its critical value. The function H3 is 
similar to Hi, and the function H4 to H2, except that H4 decreases monotonically with increasing x for all values of 
9. At the same time the range where the gauge field functions differ considerably from zero, decreases with increasing 
a and vanishes as a tends to acr- However, in this limit the gauge field functions are not continuous at the origin. 
Thus, convergence to the gauge field functions of the embedded RN solution. Hi = 0, is again pointwise. 

The Higgs field functions $1 and $2 are shown in Figs. 9h-i. <I>i increases monotonically with increasing x. ^2 is 
negative and possesses a minimum, whose position decreases with increasing a and tends to zero, as a tends to acr. 
Again the height of the extremum depends only weakly on a. Since for an embedded RN solution $1 = 1 and <I>2 — 0, 
we observe that the Higgs field functions <I>i and $2 deviate from their respective RN values in a decreasing range, as 
a approaches acr. Again, in this limit, the Higgs field functions are not continuous at the origin, thus convergence to 
the Higgs field function of the embedded RN solution is also pointwise. 



B. Coordinate transformation 



Let us focus now on the limit a acr, where the non-abelian multimonopole solutions approach the limiting RN 
solution. We recall that for n = 1 in Schwarzschild-likc coordinates a degenerate horizon forms at in = ctcr, as a 
tends to acr- For i > in the limiting solution is identical to the embedded RN solution, whereas for i < in the 
limiting solution retains its non-abelian features, and differs from any embedded abelian solution. Analogously, for 
solutions with magnetic charge n one expects that, in Schwarzschild-like coordinates, a degenerate horizon should 
form at in = riacr as a tends to acr, and that the limiting solution should retain its non-abelian features for i < in- 

In isotropic coordinates, however, the horizon radius of an extremal RN black hole is given by xn = 0. Thus, the 
limiting solution coincides in the limit a — s- acr with the extremal RN solution on the whole interval < a; < 00. 
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Therefore the question arises, as to what happens to the region, where the hmiting solution is essentially non-abelian, 
since this region shrinks to zero size in isotropic coordinates in the limit a acr- 

To elucidate this point we introduce the auxiliary Schwarzschild-like coordinate x = x^mj j , which coincides with 
the Schwarzschild-like coordinate x only, if the metric functions are independent of d. As observed above, in the 
limit a — > ctcr, the function /(a;) tends to zero on an interval of the coordinate a;, whose length tends also to zero. 
Considered as a function of x however, the function /(x) tends to zero on an interval of almost constant length. In 
this sense, the auxiliary coordinate x serves as a Schwarzschild-like coordinate. 

In Figs. lOa-c we present the n = 2 solutions in the BPS limit for a = 1.493, 1.496 and 1.499 on the interval 
< i < 4, extending beyond the critical value xh = 2acr ~ 3. As a tends to acr, the metric function / tends to 
zero on the interval < i < 3, and the metric functions m and I tend to nontrivial limiting functions. For i: > 3 the 
metric functions assume the form of the extremal RN solution with charge P = 2 and horizon radius xn — 2acr for 
a Q!cr- The behaviour of the gauge field functions and the Higgs functions is similar. For i > 3 they tend to the 
functions of the extremal RN solutions, whereas for < a; < 3 they tend to nontrivial limiting functions. 

We thus conclude, that there exists a limiting solution for a —^ aa-. This limiting solution is non-trivial and 
angle-dependent on the interior interval < i < 3. In Figs. 9b-c this limiting solution has already been reached 
in the innermost part of the interval < x < 3. For a; > 3, the limiting solution is an extremal RN solution. 
Therefore the functions of the limiting solution, and in particular the metric functions, are spherically symmetric. 
Hence, we can identify the coordinate x with the Schwarzschild-like coordinate x for i > 3, and we conclude that in 
Schwarzschild-like coordinates the horizon radius is xu = a^P, in accordance with our expectation. 

C. Gravitationally bound monopoles 

Let us now consider the mass per unit charge of the (multi)monopole solutions, to show, that gravitationally bound 
monopoles exist. The mass per unit charge of the (multi)monopole solutions decreases with increasing a and merges 
with the mass of the RN solution at acvin). In the BPS limit, for a — Q the mass per unit charge of the (multi)monopole 
solutions is precisely equal to the mass of the n = 1 monopole. For a > 0, however, we observe that the mass per unit 
charge of the multimonopoles is smaller than the mass of the n — 1 monopole. In particular, the mass per unit charge 
decreases with increasing n. Thus in the BPS limit, there is an attractive phase between like monopoles, not present 
in flat space. Moreover, multimonopoles exist for values of the gravitational coupling strength, too large for n = I 
monopoles to exist, since acrin) increases with n. The mass of the n = I monopole and the mass per unit charge of 
n = 2 and n = 3 multimonopoles in the BPS limit are shown in Fig. 11. 

For finite Higgs selfcoupling, the flat space multimonopoles have higher mass per unit charge than the n — 1 
monopole, allowing only for a repulsive phase between like monopoles. By continuity, this repulsive phase persists in 
the presence of gravity for small values of a, but it can give way to an attractive phase for larger values of a. Thus 
the repulsion between like monopoles can be overcome for small Higgs selfcoupling by sufficiently strong gravitational 
attraction. 

To mark the region in parameter space, where an attractive phase exists, we introduce the equilibrium value 
Q!oq("-i = "-2)1 where the mass per unit charge of the charge ni solution and the mass per unit charge of the charge 
n2 solution equal one another. Since amax ~ Q^cr decreases with increasing Higgs selfcoupling, and acq increases 
with increasing Higgs selfcoupling, the region of parameter space, where an attractive phase exists, decreases with 
increasing Higgs selfcoupling. In particular, for large Higgs selfcoupling, only a repulsive phase is left. 

We show the equilibrium values acq in Fig. 12. Besides the values acq(l — 2) and acq(l — 3), for which the 
monopole mass and the mass per unit charge of the n = 2 and n — S multimonopole equal one another, respectively, 
Fig. 12 also shows the values of acq(2 = 3), where the mass per unit charge of the n = 2 and the mass per unit charge 
of the n = 3 multimonopole equal one another. Thus Fig. 12 exhibits the small domain of the a-/3-plane where an 
attractive phase for like monopoles exists. 

While n = 1 monopole solutions are stable, stability of the static axially symmetric multimonopole solutions is 
not obvious. We conjecture, that the n = 2 multimonopole solutions are stable, as long as their mass per unit 
charge is lower than the mass of the n — 1 monopole. For topological number n > 3, however, solutions with only 
discrete symmetry exist in flat space j^], which, by continuity, should also be present in curved space (at least for 
small gravitational strength). For a given topological number n > 2, such multimonopole solutions without rotational 
symmetry may possess a lower mass than the corresponding axially symmetric solutions. The axially symmetric 
solutions may therefore not represent global minima in their respective topological sectors, even if their mass per unit 
charge is lower than the mass of the n = 1 monopole. 
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V. BLACK HOLE SOLUTIONS 



Here we present the static axially symmetric hairy black hole solutions of EYMH theory. We describe their 
properties, starting with the structure of the energy density of the matter fields and the deformation of the regular 
horizon. We then discuss the domain of existence of the hairy black hole solutions and describe the convergence of 
these solutions to RN solutions in the limit xa ^ x/\_ci- 

The hairy black hole solutions are then considered with respect to the results of the isolated horizon formalism. In 
particular, the mass formula is verified, and the bound state interpretation is investigated. Finally the "quasilocal 
uniqueness conjecture" is addressed. 

For the static axially symmetric hairy black hole solutions of EYMH theory we employ the same numerical technique 
as for the globally regular multimonopolc solutions (see Appendix B). The black hole solutions depend on the horizon 
radius xn and on the coupling constants a and f3. Here we consider only the BPS limit, {3 = 0. 



A. Energy density and horizon 

Let us begin the discussion of the static axially symmetric black hole solutions by considering the energy density 
of the matter fields, e. The energy density has a pronounced angle-dependence with a maximum on the p-axis. In 
particular, the energy density is not constant at the horizon. Let us now consider two representative examples for the 
energy density of the matter fields e. 

In Figs. 13a-d we exhibit the energy density of the matter fields of the n — 2 black hole solution with area parameter 
XA — 1 for a — I. Fig. 13a shows a 3-dimensional plot of the energy density as a function of the coordinates p = xsinO 
and z — xcos9 together with a contour plot, and Figs. 13b-d show surfaces of constant energy density. For small 
values of e, the surfaces of constant energy density appear ellipsoidal, being flatter at the poles than in the equatorial 
plane. With increasing values of e a toruslike shape appears. 

For smaller values of the parameters a; a and a, the energy density of the matter fields exhibits a more complicated 
structure, as seen in Figs. 14a-e, where we show the energy density of the matter fields of the n = 3 black hole solution, 
with area parameter xa — 0.5 for a — 0.5. Fig. 14a again shows a 3-dimensional plot of the energy density as a 
function of the coordinates p = a;sin6' and z = xcos9 together with a contour plot, while Figs. 14b-e show surfaces 
of constant energy density. Whereas the surfaces of constant energy density still appear ellipsoidal for small values 
of e, here with increasing values of e the toruslike shape appears together with two additional ellipsoids covering the 
poles. For the largest values of the energy density only the toruslikc shape remains. 

The n-dependence of the energy density of the matter fields is illustrated in Fig. 15, where we show the energy 
density of the black hole solutions with magnetic charges n = 1 — 3 and area parameter xa = 0.5 for a ~ 0.5. With 
increasing magnetic charge n the absolute maximum of the energy density of the solutions, residing on the p-axis, 
shifts outward and decreases significantly in height. 

Let us now turn to the regular horizon of the hairy black hole solutions, which resides at a surface of constant 
radial coordinate x = xjj. Even though the radial coordinate is constant at the horizon, the horizon is deformed. The 
deformation is revealed, when measuring the circumference of the horizon along the equator, Le, and the circumference 
of the horizon along the poles, Lp, 




Lp — 2 



x—x^,0—7t/2 



d9 . / — X 



(75) 



x—xnj^p— const. 



since the hairy black hole solutions have Lp ^ (in general). The deviation from spherical symmetry is small, 
though. For instance, for the solution of Fig. 13 L^/Lp ~ 0.99076, and for the solution of Fig. 14 L^/Lp — 0.99977. 

The hairy black holes satisfy the zeroth law of black hole mechanics, which states that the surface gravity n is 
constant on the horizon |^^. This is dictated by the full set of EYMH equations, as discussed in section II. D. 
Numerically the surface gravity is also constant, as demonstrated in Fig. 16 for the n = 2 solution with horizon area 
parameter xa = 1 for a = 1. In constrast to the surface gravity itself, the expansion coefficients f2{d) and m2{0) of 
the metric functions / and to, which enter into the expression for k, Eq. (^2|), possess a non-trivial angular dependence 
at the horizon, as seen in Fig. 16. 
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B. Domain of existence 



The domain of existence of the axially symmetric non-abeUan black holes is very similar to the domain of existence 
of the spherically symmetric black hole solutions. For a fixed value of the coupling constant a, a < Q!max(^), hairy 
black hole solutions emerge from the globally regular solution in the limit xa 0. They persist up to a maximal 
value of the area parameter a;A,max, which depends on n and limits the domain of existence of hairy black holes. The 
domain of existence of hairy black hole solutions with magnetic charge n = 2 is indicated in Fig. 3, where crosses 
mark the maximal values XA,max obtained. 

The domain of existence of the hairy black hole solutions with n > 1 also consists of two regions with distinct 
critical behaviour, separated by a particular value of a, denoted a{n). For a{n) < a < amaxin), the hairy black hole 
solutions bifurcate at a critical value XA.cr with an extremal RN solution with charge n, while for < a < a{n), the 
hairy black hole solutions bifurcate at a critical value XA.cr with a non-extremal RN solution with charge n. For small 
values of a, the branch of hairy black hole solutions extends backwards from the maximal value a;A,max to the critical 
value a;A,cr, whereas for larger values of a, both values coincide, XA,max = xa,ct- This pattern of the n > 1 hairy 
black hole solutions is completely analogously to the pattern observed for the n = 1 solutions. Our numerical results, 
as exhibited in Fig. 3, are consistent with the conjecture d(2) = •\/3/2, suggesting that a{n) is independent of n. 

C. Area parameter a^A-dependence 

Before demonstrating the dependence of the hairy black hole solutions on the area parameter .ta and the convergence 
of the solutions to the limiting RN solution, we briefly consider the relation between the area parameter xa and the 
horizon radius in isotropic coordinates xh, since Xn is one of the parameters employed in the calculations. 

In the limit Xu ^ 0, hairy black hole solutions emerge from the globally regular solution. With increasing parameter 
xh the horizon area of the hairy black hole solutions increases, until a maximal value of the parameter xu is reached. 
This maximal value of the parameter xu, however, does not coincide with the maximal value of the horizon area, and 
thus the maximal value of the area parameter xa- In the further discussion we consider the two regions of the domain 
of existence with distinct critical behaviour separately. 

For small values of a, < a < ain), the branch of hairy black hole solutions, as a fimction of the horizon radius 
Xu, extends backwards from the maximal value a;H,max, until the critical value xn.cr of the bifurcation with the non- 
extremal RN solution is reached. When xu decreases from xn.max, the area parameter increases further up to its 
maximal value .TA.max, from where it decreases, until it reaches the critical value xa.ct of the bifurcation. This is 
illustrated in Fig. 17 for the n = 2 black hole solutions for a = 0.5. The endpoint of the curve marks the bifurcation 
with the non-extremal RN solution at a finite value of the area. 

For larger values of a, a{n) < a < ama^in), the branch of hairy black hole solutions bifurcates with an extremal 
RN solution, which has Xu = 0. As a function of the horizon radius xu, the branch of hairy black hole solutions also 
extc;n(ls backwards from the maximal value a;H,max, but it ends at the critical value xu,ci = 0, where it bifurcates with 
the extremal RN solution. The area parameter, in contrast, increases monotonically, and reaches its maximal value 
a^A.max = XA,cr Sit the bifurcation, as illustrated in Fig. 17 for the n = 2 black hole solutions for a = 1. 

Let us now consider the dependence of the hairy black hole solutions on the area parameter xa in more detail, 
distinguishing again the regions < a < d(n) and a{n) < a < amax('^)- For a{n) < a < amaxin), the limiting 
solution of the branch of hairy black hole solutions is an extremal RN solution. Since convergence to an extremal 
RN solution has been considered in great detail in section IV, in the following we put more emphasis on the region 
< a < d(n), where convergence to a non-extremal RN solution is observed. 

For small values of a in the region < a < a(n), the branch of hairy black hole solutions extends backwards from 
the maximal value a;A,max, which limits the domain of existence, to the critical value XA,cr, where the bifurcation 
occurs. In Figs.l8a-c we exhibit a set of n = 2 black hole solutions for four values of the area parameter along the 
black hole branch. Increasing xa from zero along the branch of hairy black hole solutions, first the solutions with 
Xa = 0.46 and 1.09 are encountered, then the solution with xa = 1-39 close to the maximal value a;A,max is passed, 
and finally the solution with xa = 1-09 close to the critical value xa,ci is reached. 

The energy density of the matter fields is shown in Fig. 18a, the metric function / in Fig. 18b, and the norm of 
the Higgs field | $ |= y^f + ^f in Fig. 18c. All functions change drastically along the hairy black hole branch. In 
particularly, we observe, that they approach the limiting functions of the corresponding non-extremal RN solution for 
Xa — a^A,cr ~ 1.07. Thus in this limit the functions become spherically symmetric. 

Convergence to the limiting RN solution is nicely seen in Fig. 19, where H2{xii) and ^i{xii) are shown, the horizon 
values of the gauge field function H2 and the Higgs field function $1, respectively. Starting from the unique values 
1^2(0) = 1 and $1(0) = of the regular solution at the origin, for finite xa an angle-dependence arises, which becomes 
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maximal as the maximal value of the area parameter a; a, max is approached. On the backward bending part of the 
hairy black hole branch, the angle-dependence diminishes and vanishes in the limit xa — a;A,cr, where the horizon 
values tend to the unique horizon values of the corresponding RN solution, H2{x}i) — and <l>i(a;H) — 1- 

For a > d, the limiting solution is reached at the maximal value of the area parameter a;A,max = 2:A,cr- Since 
a^A.cr < na, a gap between the branch of hairy black hole solutions and the corresponding RN branch occurs, 
similarly as in the case of the globally regular solutions. 

Let us now consider the mass, the surface gravity and the deformation of the horizon of the hairy black hole 
solutions. In Fig. 20 we see the dependence of the mass of the n = 2 black hole solutions on the area parameter 
xa for a = 0.5 and a = 1. (Fig. 20 is completely analogous to Fig. 4, where the mass of the n = 1 black hole 
solutions is shown.) For a = 0.5 < a, the mass increases with increasing xa, reaching its maximal value at a;A.max- 
Forming a spike, the mass then decreases with decreasing xa, and reaches the mass of the limiting non-extremal 
RN solution with charge n = 2 at XA,cr- For comparison, the corresponding branch of RN solutions is also shown. 
Apart from the critical point XA,cr, the RN branch and the hairy black hole branch also coincide at the crossing point 
2:a, cross ~ 1.32. Consequently, the RN solutions have a lower mass than the non-abelian black hole solutions, when 
2^A, cross < XA < XA,ci- For a = 1 > a, the mass increases monotonically with increasing xa, reaching its maximal 
value at the maximal value of the area parameter XA.max = a^A,cr- Since XA.cr < 2a, a gap between the branch of 
hairy black hole solutions and the corresponding RN branch occurs. In Fig. 20, the RN branch would begin at xa = 2. 

In Fig. 21 the dependence of the surface gravity of the n = 2 hairy black hole solutions on the area parameter xa 
is shown for a = 0.5 and a — 1. (Fig. 21 is completely analogous to Fig. 5, where the surface gravity of the n — 1 
black hole solutions is shown.) Again, convergence the corresponding RN values is observed. 

Let us now turn to the shape of the horizon. In Fig. 22 the dependence of the ratio of circumferences L^/Lp for the 
n = 2 black hole solutions on the area parameter xa is shown for a = 0.5 and a = 1. Interestingly, for a — 0.5 < a 
the ratio L^/Lp, shows a rather complicated dependence on the area parameter. With increasing area parameter 
the ratio L^/Lp first decreases, until it reaches a minimum close to XA.max- Along this part of the branch of hairy 
black hole solutions Le/Lp < 1. Continuing along the hairy black hole branch, the ratio then increases, first with 
increasing xa until XA.max is reached, and then with decreasing xa- Passing Le/Lp = 1, the ratio increases further 
with decreasing a; a, and attains a maximum L^/Lp > 1. From there L^/Lp decreases with decreasing xa and reaches 
the limiting value L^/Lp = 1 of the spherically symmetric RN solution for xa — > a^A.cr- Along the latter part of 
the branch of hairy black hole solutions, L^/Lp > 1. The occurrence of L^/Lp > 1 is a new phenomenon, not seen 
previously. (The EYM solutions had only values of the ratio L^/Lp < 1 [TtI].) 

For a = \ > a, the dependence of the ratio L^jLp on the area parameter is much simpler. Again, with increasing 
area parameter the ratio L^jLp decreases, until it reaches a minimum not too far from XA.max- The ratio then 
increases with increasing xa, to reach the limiting value L^/Lp = 1 of the spherically symmetric RN solution for 
XA XA,cr- Thus for a — 1 the hairy black hole solutions always have Le/Lp < 1. We also note that the magnitude 
of the deformation of the horizon of the a — 1 solutions is considerably greater than the magnitude of the deformation 
of the horizon of the a — 0.5 solutions. 

Let us finally consider the dependence of the deformation of the horizon on the magnetic charge n. For the solutions 
with charge n > 2 we observe the same features of the ratio Le/Lp as for the solutions with charge n = 2, described 
above. This is seen in Fig. 22, where the dependence of the ratio of circumferences Le/Lp for the n = 3 black hole 
solutions on the area parameter xa for a = 0.5 < d is also shown. We observe, that the maximal deformation of the 
horizon increases with n. 



D. Isolated horizon results 

Let us now address the predictions of the isolated horizon formalism for the axially symmetric hairy black hole 
solutions with n > 1. We first consider relation (|5j) between black hole mass, soliton mass and horizon mass for 
the n > 1 solutions. Our numerical calculations confirm this relation for the axially symmetric hairy black hole 
solutions of EYMH theory. In Fig. 23 we show the dependence of the horizon mass per unit charge ^jla/ {o^n), 
obtained according to Eq. (^5|), on the area parameter xa for the n — 2 black hole solutions for a = 0.5 and a = 1. 
Adding the mass per unit charge of the corresponding multimonopole solutions (/ircg/(2Q;^)(Q; — 0.5) — 0.961105 and 
/ircg/ (2q!^)(q! = 1.0) — 0.847943, respectively,) precisely gives the values of the mass per unit charge of the hairy black 
hole solutions in Fig. 20. 

Next we consider the interpretation of the hairy black holes as a bound states of regular solutions and Schwarzschild 
black holes |Q. We therefore identify the binding energy of these systems /Xbind/a^, according to Eq. (|57|). In Fig. 24 
we present the dependence of the binding energy /ibind/a^ on the area parameter xa for the n = 2 black hole solutions 
for a = 0.5 and a = 1. For comparison, the binding energy of the corresponding RN solutions is also given. Fig. 24 
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is completely analogous to Fig. 7, where the binding energy of the n = 1 black hole solutions is shown. 

However, for the n > 1 hairy black hole solutions, we can also consider other compound systems with magnetic 
charge n, consisting of a soliton and a black hole. Therefore the notion of binding energy is no longer unique, and 
the binding energy of the compound system of a, n — 2 multimonopole and a Schwarzschild black hole represents only 
one particular case. 

In the following we restrict the discussion to hairy black hole solutions with magnetic charge n — 2. Let us define 
the binding energy /ibind/'^^ ^^'^ compound system of a n = 2 multimonopole and a Schwarzschild black hole by 

and compare it to the binding energy Mbind/"^^ ^he compound system of two n = 1 monopoles and a Schwarzschild 
black hole, 

;,(»=2) (a,^) = 2^(^=1) + ^,sixA) + mLcI - + MLd ■ (77) 

Since the binding energy /ibind/*^^ differs from the binding energy fi^^-^^^/a^ by the difference between the mass of the 
n = 2 multimonopole and the mass of two n — 1 monopoles, the binding energy fj.'^^^^/a^ is smaller than the binding 
energy Mbind/*^^' shown in Fig. 24. 

Alternatively we can consider a n = 2 hairy black hole as a compound system of a n = 1 monopole and a charged 
black hole. For the compound system of a n = 1 monopole and a ?i = 1 hairy black hole, the binding energy Aibind/*^^ 
is given by 

/i^"='H^A) = + ^^^''='\XA) + M^ind = + A^bind , (78) 

in the parameter range, where both black hole solutions coexist. Likewise, for the compound system of a n = 1 
monopole and a RN black hole with unit charge, the binding energy Mbind/*^^ given by 

/i("=2) {xa) = M^og"^^ + MRn(2:a) + Aibind = /^'* + A^bind , (79) 

in the parameter range, where both solutions coexist. 

Let us now illustrate these possibilities for the binding energy by considering the masses of the corresponding 
compound systems. In Figs. 25a-b we show the dependence the masses /i^/a^ — /i^/a^, Eqs. (p6|)-(^9|), on the area 
parameter xa, and compare it to the mass of the n = 2 hairy black hole solutions, ^jS^'^'^^x a) / cP' , for a = 0.5 and 
a = 1, respectively. 

The binding energy jjb^-^^^/a^ and the binding energy ^bind/*^^ ^'^s always negative for the n = 2 hairy black hole 
solutions. In contrast, the binding energy Mbi„(j/a^, representing the binding energy of the compound system of a 
n = 1 monopole and a ti = 1 hairy black hole, changes sign. For a = 0.5, Mbind/*^^ positive beyond xa ~ 0.34, 
and for a = 1, it is positive beyond xa ~ 0.52. The binding energy Mbind/'^^' representing the binding energy of the 
compound system of a n = 1 monopole and a RN black hole with unit charge, changes sign at xa ~ 0.63 for a = 0.5, 
whereas it is always positive for a = 1. 

Thus for a = 0.5, the n = 2 hairy black holes possess a mass lower than the masses /i^ — /i* of the compound 
systems only for xa < 0.34. For xa > 0.34, the n — 2 hairy black hole is either heavier than the compound system 
of a n = 1 monopole and a n = 1 hairy black hole, or (when the n — I hairy black hole branch ceases to exist) it is 
heavier than the compound system of a rt = 1 monopole and a RN black hole with unit charge. 

For a — 1 the n ^ 2 hairy black hole is energetically favourable for xa < 0.52, where the compound system of a 
n — 1 monopole and a rt = 1 hairy black hole is heavier, and for XA,cr(n = 1) < < 1, since below xa = 1 no RN 
solutions exist. {xA.ci{n = 1) < xa < 1 represents the gap between the n = 1 hairy black hole branch and the RN 
branch with unit charge.) 

Let us finally consider the "quasilocal uniqueness conjecture" , which states that static black holes are uniquely 
specified by their horizon area and horizon charges. In Fig. 26 we show the dependence of the non-abelian horizon 
magnetic charge on the area parameter xa for the n — 2 hairy black hole solutions for a = 0.5 and 1. The horizon 
charge increases monotonically along the branches of n = 2 hairy black hole solutions, analogously to the horizon 
charge of the n ~ 1 hairy black hole solutions, shown in Fig. 8. In the limit xa ^ a^A.cr, for a = 0.5 the value of 
the RN solution with charge n = 2 is reached, as expected, whereas for a — 1 the limiting value is smaller than two 
(reflecting the gap between the hairy black hole branch and the RN branch). The horizon electric charge vanishes. 

To address the "quasilocal uniqueness conjecture", we now also consider the branches of embedded RN solutions 
with integer values of the magnetic charge, beginning at xa = ctn. The lowest branch of embedded RN solutions has 
unit magnetic charge, and thus unit horizon magnetic charge. For a = 0.5, this branch begins at xa = 0.5, and crosses 
the corresponding (a = 0.5) branch of rt = 2 hairy black hole solutions. Consequently, at the crossing point there 
exist two distinct black hole solutions, one with hair and the other without hair, with the same horizon area and the 
same non-abelian horizon magnetic charge, representing a counterexample to the "quasilocal uniqueness conjecture". 
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VI. CONCLUSIONS 



We have considered static axially symmetric multimonopole and black hole solutions in EYMH theory. We have 
presented these solutions in detail and discussed their properties. Our particular interests were the investigation of 
an attractive phase between like monopoles and the study of monopole and black hole properties predicted by the 
isolated horizon formalism. 

Concerning the multimonopole solutions we observe, that the mass per unit charge of the (multi) monopole solutions 
decreases with increasing a. In the BPS limit, for a = the mass per unit charge is precisely equal to the mass of the 
n = 1 monopole. For a > 0, however, we observe that the mass per unit charge of the multimonopoles is smaller than 
the mass of the n = 1 monopole. In particular, the mass per unit charge decreases with increasing n. Thus in the 
BPS limit, there is an attractive phase between like monopoles, not present in flat space. Moreover, multimonopoles 
exist for gravitational coupling strength, too large iov n — 1 monopoles to exist. 

For finite Higgs selfcoupling, the flat space multimonopoles have higher mass per unit charge than the n — 1 
monopole, allowing only for a repulsive phase between like monopoles. By continuity, this repulsive phase persists in 
the presence of gravity for small values of a, but it can give way to an attractive phase for larger values of a. Thus 
the repulsion between like monopoles can be overcome for small Higgs selfcoupling by sufficiently strong gravitational 
attraction. At the equilibrium value aeq, multimonopole mass per unit charge und monopole mass equal one another. 
The equilibrium value ckcq increases with increasing Higgs selfcoupling, yielding a decreasing region in parameter space 
for the attractive phase. For large Higgs selfcoupling, only a repulsive phase is left. 

While singly charged monopole solutions are stable, stability of the static axially symmetric multimonopole solutions 
is not obvious. We conjecture, that the n = 2 multimonopole solutions are stable, as long as their mass per unit 
charge is lower than the mass of the n = 1 monopole. For topological number n > 3, however, solutions with only 
discrete symmetry exist in flat space ||], which, by continuity, should also be present in curved space (at least for 
small gravitational strength). For a given topological number n > 2, such multimonopole solutions without rotational 
symmetry may possess a lower mass than the corresponding axially symmetric solutions, as suggested by analogy from 
the multiskyrmions in flat space . The axially symmetric solutions may therefore not represent global minima in 
their respective topological sectors, even if their mass per unit charge is lower than the mass of the n = 1 monopole. 

Let us now turn to the black hole solutions of SU(2) EYMH theory. Besides embedded abelian black hole solutions, 
SU(2) EYMH theory also possesses genuine non-abelian black hole solutions |l^. The static SU(2) EYMH black hole 
solutions are no longer uniquely determined by their mass and charge alone. Indeed, in a certain region of the domain 
of existence of hairy black hole solutions, also embedded RN solutions with the same mass and charge exist. The 
non-abelian black hole solutions therefore represent counterexamples to the "no- hair" conjecture jT^ . 

While static spherically symmetric {n = 1) EYMH black holes were studied in great detail non-abelian black 
hole solutions with magnetic charge n > 1 were previously only considered perturbatively jl^ . We have obtained static 
axially symmetric black hole solutions with integer magnetic charge n > 1 numerically. These black hole solutions 
are asymptotically flat, and they possess a regular deformed horizon. Being static and not spherically symmetric, 
these black hole solutions represent further examples, showing that Israel's theorem cannot be generalized to EYM 
or EYMH theory. While previous (non-perturbative) counterexamples ||2^,^ were classically unstable, hairy EYMH 
black holes should provide classically stable counterexamples . 

Considering the static axially symmetric solutions from the isolated horizon formalism point of view, we have verified 
the mass relation between the monopole and the black hole solutions, showing that the black hole mass is given by 
the sum of the soliton mass and the horizon mass. Interpreting the hairy black holes as bound states of solitons and 
Schwarzschild black holes ||2^, we have studied the binding energy of these bound systems. We have furthermore 
considered the binding energy with respect to various other compound systems, such as a n — 1 soliton and a n = 1 
hairy black hole or a n — 1 soliton and a RN black hole with unit charge. The "quasilocal uniqueness conjecture" 
claims, that black holes are uniquely specified by their horizon area and their horizon electric and magnetic charge. 
Since we have constructed a counterexample to this conjecture, the need for a new formulation of the uniqueness 
conjecture arises. 

The hairy black hole solutions studied here represent only the simplest type of non-spherical black hole solutions. 
Indeed, there are gravitating black hole solutions with much more complex shapes and only discrete symmetries left. 
In curved space such black holes without rotational symmetry, have been considered so far only perturbatively . It 
remains a challenge to construct such solutions non-perturbatively and to find out, whether such black hole solutions 
without rotational symmetry are stable. 
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VII. APPENDIX A 



A. Tensors F^„, D^<^, 



We expand the field strength tensor and the covariant derivative with respect to the PauU matrices 
^, see Eq. ® ) 



and 



Inserting ansatz ( |l^ for the gauge field, we obtain the non-vanishing coefficients Fyi^-* and D^^'^-^^ 



F^::^^ = --{H,,g + rH2,r) , 
r 



sm 



= n +H1H3 + cote Hi) 



i^i;^ = -n sin 9 {H3,e - 1 + H2H4 + cotOH^) 



F^^J = nsine (i/4,e - H2H3 - cote {H2 - Hi)) , 



and F^^) = 



1 



1 , 

= - r$2,r) , 

r 

= ($2,9 +i?2*l) , 

= nsin^ (iJ4$i + i?3$2 + cot6'<i>2) 



It is convenient to define 



^2 

rip 



Fl 



{rHi^r~H2,e) , 



where the second term in the definition of F^g represents the gauge fixing term, as well as 
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With the ansatz for the metric (|l^) we obtain the Lagrange densities 

1 



Lf 

i<[> 



2m I r^m ^ sin^ 0/ ^ ^2 ^e^- 



2 m 



1 



and the non-vanishing components of the stress energy tensor, 



(84) 
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sin Ql 



2 



1 



r'' ir^ sm 



m A 
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.2 ^^^-^a 



Ir^ sm 



/ 8 
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/ 8- 



i fi^2 

sm ^ 



Zr2sin^ 61 



Zr2 sin2 61 A 
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VIII. APPENDIX B 



Subject to the corresponding set of boundary conditions, we solve the system of coupled non-linear partial differential 
equations numerically. For the globally regular solutions, we employ the radial coordinate 

(86) 



instead of x, to map the infinite interval of the variable x onto the finite interval [0, 1] of the variable z. For the 
derivatives this leads to the substitutions 

rF^r — > z{l - z)F^ , 
r''F^r,r ((1 - z)2f,,, - 2(1 - z)F,) (87) 

(88) 

for any function F in the differential equations. In this form we have solved the system of differential equations 
numerically. 

To map spatial infinity to the finite value z = 1, we employ for the black hole solutions the radial coordinate 

z - 1 - — . (89) 

X 

For the derivatives this leads to the substitutions 

ri^,, (1 - z)F, , 
r^F^r,r — > (1 - zfF_,_,,_ - 2(1 - z)F^ (90) 

(91) 

for any function F in the differential equations. 

For the black hole solutions we furthermore introduce the functions /(z, 0), ?7i(z, 0) and l{z, 6) ||l^ , [l7|] , 

/-(z, ^) = ^ , ^(z, 0)^^, /-(z, 0) = (92) 
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where z is the compactified coordinate (p9|). Since in the hmit x ^ cx), the variable z approaches the value 1, the 
boundary conditions for the functions /, to and I coincide with the boundary conditions for the functions /, to and I 
at infinity. At the horizon, the boundary conditions of the functions /, to and I are given by 

(/-aJ)U=o = 0, (m + 9,TO)U=o =0 , (/"+aJ)|,^o = . (93) 

To satisfy the regularity condition (|ll|) in the numerical calculations, we have introduced the new function 17(2, d) 

9i^^0)^^ (94) 

with the boundary conditions on the symmetry axis and at the horizon 

9\e=o = 1 , d,g\,=o = . (95) 

The numerical calculations are performed with help of the program FIDISOL, which is extensively documented 
in iQ. The equations are discretized on a non-equidistant grid in z and 9. Typical grids used have sizes 150 x 30, 
covering the integration region < 2; < 1 and < < tt/2. 

The numerical method is based on the Newton- Raphson method, an iterative procedure to find a good approximation 
to the exact solution. Let us put the partial differential equations into the form P{u) = 0, where u denotes the unknown 
functions (and their derivatives). For an approximate solution u^^\ P{u^^'>) does not vanish. If we could find a small 
correction Au, such that u^^^ + Au is the exact solution, P{u^^'> + Am) = should hold. Approximately the expansion 
in Au gives 

BP 

= P(u(i) + Au) » P(u(i)) + ^(u^^hAu . 

au 

The equation P{u'^^^) = — ^(u^^-')Au can be solved to determine the correction Au^^^ = Au. u^^^ = u^^^ + Au'^-'^^ 
will not be the exact solution but an improved approximate solution. Repeating the calculations iteratively, the 
approximate solutions will converge to the exact solutuion, provided the initial guess solution is close enough to the 
exact solution. The iteration stops after i steps if the Newton residual P(u('^) is smaller than a prescribed tolerance. 
Therefore it is essential to have a good first guess, to start the iteration procedure. Our strategy therefore is to use a 
known solution as guess and then vary some parameter to produce the next solution. 

To construct axially symmetric EYMH solutions, we have used the known spherically symmetric EYMH solutions 
as starting solutions with n = 1. We have then increased the 'parameter' n slowly, to obtain the desired axially 
symmetric solutions at integer values of n. 

For a numerical solution it is important to have information about its quality, i. e. to have an error estimate. The 
error originates from the discretization of the system of partial differential equations. It depends on the number 
of gridpoints and on the order of consistency of the differential formulae for the derivatives. FIDISOL provides an 
error estimate for each unknown function, corresponding to the maximum of the discretization error divided by the 
maximum of the function. For the solutions presented here the estimations of the relative error for the functions are 
on the order of lO^'^ for n — 2 and n = 3. 
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metric function N 




FIG. la. The metric function A' of the n = 1 monopolc solution is shown as a function of the compactified dimensionless 
Schwarzschild-likc coordinate Zs = x/{l + x) in the BPS limit for five values of a along the monopole branch, in particular for 
a value close to the maximal value of a, amax ~ 1.403 and for a value close to the critical value of a, acr « 1.386. 
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FIG. lb. The same as Fig. la for the gauge field function K. 
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FIG. Ic. The same as Fig. la for the Higgs field function H. 



metric function f 

1.0 
0.8 

0.6 
0.4 

0.2 

0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

z 

FIG. 2a. The metric function / of the n = 1 monopole solution is shown as a function of the compactifled dimensionless 
isotropic coordinate z = x / (1 + x) m the BPS limit for five values of a along the monopole branch, in particular for a value 
close to the maximal value of a, Omax ~ 1.403 and for a value close to the critical value of a, Ocr ~ 1.386. 
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FIG. 2b. The same as Fig. 2a for the gauge field function K. 
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FIG. 2c. The same as Fig. 2a for the Higgs field function R. 
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FIG. 3. The domain of existence of the hairy black hole solutions in the BPS limit is shown in the i-A-Q-planc. The solid 
line shows the maximal values XA.max obtained for the n — 1 hairy black hole solutions, while the crosses represent the maximal 
values XA.max obtained for n = 2 hairy black hole solutions. The asterisk maxks the value a(2) = V3/2, conjectured to separate 
the two regions with distinct critical behaviour. Also shown are the extremal RN solutions with unit charge and charge two. 
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FIG. 4. The dependence of the mass ^j./o? of the n — 1 hairy black hole solutions on the area parameter xa is shown in 
the BPS limit for a = 0.5 and a = 1. For comparison, the mass of the corresponding RN solutions is also shown. 
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FIG. 5. The dependence of the surface gravity k of the n = 1 hairy black hole solutions on the area parameter xa is shown 
in the BPS limit for a = 0.5 and a = 1. For comparison, the surface gravity of the corresponding RN solutions is also shown. 
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FIG. 6. The dependence of the horizon mass ha. jo? of the n = 1 hairy black hole solutions on the area parameter xa is 
shown in the BPS limit for a = 0.5 and a = 1. 
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FIG. 7. The dependence of the binding energy /ibind/a^ of the n = 1 hairy black hole solutions on the area parameter xa 
is shown in the BPS limit for a = 0.5 and a = 1. For comparison, the binding energy of the corresponding RN solutions is also 
shown. 
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FIG. 8. The dependence of the non-abelian horizon magnetic charge P{xa) of the n = 1 hairy black hole solutions on the 
area parameter xa is shown in the BPS limit for a = 0.5 and a = 1. 
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FIG. 9a. The metric function / of the axially symmetric n = 2 monopole solution is shown as a function of the compactified 
dimensionless isotropic coordinate z = x/{l + x) for the angles = = ■n/A and 6 = 7r/2 in the BPS limit for four values of 
a along the multimonopole branch, a = 1, 1.2, 1.4 and 1.499. 
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FIG. 9b. The same as Fig. 9a for the metric function m. 
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FIG. 9c. The same as Fig. 9a for the metric function I. 
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FIG. 9d. The same as Fig. 9a for the gauge field function H\ . 
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FIG. 9e. The same as Fig. 9a for the gauge field function Hi- 
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FIG. 9f. The same as Fig. 9a for the gauge field function H3. 
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FIG. 9g. The same as Fig. 9a for the gauge field function H4. 
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FIG. 9h. The same as Fig. 9a for the Higgs field function 3>i . 
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FIG. 9i. The same as Fig. 9a for the Higgs field function $2- 
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FIG. 10a. The metric function / o f the axially symmetric n = 2 monopole solution is shown as a function of the auxiliary 

Schwarzschild-likc coordinate x = x ^fmj], for the angles 6 = 0, 6 = 7r/4 and Q = 7r/2 in the BPS limit for a = 1.493, 1.496 
and 1.499, close to the critical a. Also shown is the metric function / of the RN solution with charge P = 2 and area parameter 
XA = 1.5. 
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FIG. 10c. The same as Fig. 10a for the norm of the Higgs field |$| = ^/^T^. 
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FIG. 11. The dependence of the mass per unit charge ii/(a^n) of the hairy black hole solutions on the parameter a is shown 
in the BPS limit for magnetic charge n = 1, 2 and 3. For comparison, the mass per unit charge of the extremal RN solutions 
is also shown. 
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FIG. 12. The dependence of the maximal values amax(?^) of the hairy black hole solutions on the parameter (3 is shown for 

magnetic charge n = 1, 2 and 3. Also shown arc the oquihbrium values acq(ni = n2), in particular Q:cq(l = 2), Oeq(l = 3) and 
aeq(2 = 3), for which the mass per unit charge of the charge ni and and the charge n2 (multi)monopole equal one another. 
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FIG. 15. The energy density e = —Tq of the hairy black hole solutions with magnetic charge n = 1, 2 and 3 and area 
parameter xa = 0.5 is shown as a function of the dimensionless isotropic coordinate x for the angles 9 = 0, 9 = 7r/4 and 
9 = Tx/2 in the BPS limit for a = 0.5. 
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FIG. 16. The anglo-indcpondoncc of the surface gravity k is shown for the hairy black hole solution with magnetic charge 
n = 2 and area parameter .xa = 1-0 in the BPS hmit for a = 1.0. Also shown is the angle-dependence of the normalized 
expansion coefHcients /2 and m2 of the metric functions / and m. 
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FIG. 17. The dependence of the horizon radius in isotropic coordinates xn of the n = 2 hairy black hole solutions on the 
area parameter xa is shown in the BPS limit for a = 0.5 and a = 1. 
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FIG. 18a. The energy density e = — Tq of the n = 2 hairy black hole solution is shown as a function of the dimensionless 
compactified coordinate a = 1 — ^ for the angles = 0, 9 = 7r/4 and = -it/2 in the BPS limit for four values of a along the 
hairy black hole branch, in particular for a value close to the maximal value of max w 1.39, and for a value close to the 

critical value of 1.07. 
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FIG. 18b. The same as Fig. 18a for the metric function /. 



norm of Higss fr 



1.00 



^ 0.50 



0.25 . 



0.00 





1 


1 


1 1 j 

x^=1.09 / 






= 








= 7t/4 
= 7t/2 


r I \ 
/// f I 

X4=i.39.vr ./ I 

. ■ ■ //I \ 

■■■■■■y^'y/f \ 








- ^ ' 7' 








^ — A^ .•// - 


- - " ^ • ■ 'i^^ * — 




-L 


1 _ _ 





0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 18c. The same as Fig. 18a for the norm of the Higgs field | $ |= y'4>f + 
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FIG. 19. The dependence of the horizon value of the gauge field function H2{xh) and the horizon value of the Higgs field 
function ^i(xh) of the n = 2 hairy black hole solutions on the area parameter xa is shown in the BPS limit for a = 0.5. 
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FIG. 20. The dependence of the mass per unit charge /i/ (a'^n) of the n — 2 hairy black hole solutions on the area parameter 
XA is shown in the BPS limit for a = 0.5 and a = 1. For comparison, also the mass per unit charge of the RN solutions for 
a = 0.5 is shown. For a=1.0, RN solutions with charge P = 2 exist only for ka > 2. 
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FIG. 21. The dependence of the surface gravity n of the n = 2 hairy black hole solutions on the area parameter xa is shown 
in the BPS limit for a = 0.5 and a = 1. For comparison, also the surface gravity of the RN solutions for a = 0.5 is shown. 
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FIG. 22. The dependence of the ratio of circumferences Le/Lp of the n = 2 hairy black hole solutions on the area parameter 
Xa is shown in the BPS limit for a = 0.5 and a = 1. Also shown is the ratio of circumferences Le/Lp of the n = 3 hairy black 
hole solutions for a = 0.5. 
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FIG. 23. The dependence of the horizon mass ha/oc^ of the n = 2 hairy black hole solutions on the area parameter xa is 
shown in the BPS limit for a = 0.5 and a = 1. 
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FIG. 24. The dependence of the binding energy /Xbind/a^ of the n = 2 hairy black hole solutions on the area parameter xa 
is shown in the BPS limit for a = 0.5 and a = 1. For comparison, the binding energy of the corresponding RN solutions is also 
shown. 
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FIG. 25a. The dependence of the mass of the compound system of a n = 2 multimonopole and a Schwarzschild black 
hole, of two n = 1 monopolcs and a Schwarzschild black hole, of a n = 1 monopole and a n = 1 hairy black hole, and of a 
n = 1 monopole and a RN black hole with unit charge, on the area parameter xa is shown in the BPS limit for a = 0.5. For 
comparison the mass of the n = 2 hairy black hole solutions is also shown. The inlet shows the binding energy with respect to 
the compound system of a n = 1 monopole and a n = 1 hairy black hole. 




FIG. 25b. The same as Fig. 25a for a = 1. 
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FIG. 26. The dependence of the non-abehan horizon magnetic charge P{xa) of the n = 2 hairy black hole solutions on the 
area parameter xa is shown in the BPS limit for a = 0.5 and a = 1. Also shown is the horizon magnetic charge of the RN 
solutions with charge n = 1 and n = 2 for a = 0.5. 



46 



